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In this paper, we combine thermal effects with Landau-Zener (LZ) quantum tunneling effects 
in a dynamical Monte Carlo (DMC) framework to produce satisfactory magnetization curves of 
single-molecule magnet (SMM) systems. We use the giant spin approximation for SMM spins and 
consider regular lattices of SMMs with magnetic dipolar interactions (MDI). We calculate spin rever- 
sal probabilities from thermal-activated barrier hurdling, direct LZ tunneling, and thermal-assisted 
LZ tunnelings in the presence of sweeping magnetic fields. We do systematical DMC simulations for 
Mni2 systems with various temperatures and sweeping rates. Our simulations produce clear step 
structures in low-temperature magnetization curves, and our results show that the thermally acti- 
vated barrier hurdling becomes dominating at high temperature near 3K and the thermal-assisted 
tunnelings play important roles at intermediate temperature. These are consistent with corre- 
sponding experimental results on good Mni2 samples (with less disorders) in the presence of little 
misalignments between the easy axis and applied magnetic fields, and therefore our magnetization 
curves are satisfactory. Furthermore, our DMC results show that the MDI, with the thermal effects, 
have important effects on the LZ tunneling processes, but both the MDI and the LZ tunneling give 
place to the thermal-activated barrier hurdling effect in determining the magnetization curves when 
the temperature is near 3K. This DMC approach can be applicable to other SMM systems, and 
could be used to study other properties of SMM systems. 

PACS numbers: 75.75.-c, 05.10.-a, 75.78.-n, 75.10.-b, 75.90,+w 



I. INTRODUCTION 

Single-molecule magnet (SMM) systems attract more 
and more attention because they can be used to 
make devices for spintronic applications 1,2 , quan- 
tum computing 3 , high-density magnetic information 
storage^ etc£~— . Usually, a SMM can be treated 
as a large spin with strong magnetic anisotropy at 
low temperature. The most famous and typical 
is Mni2-ac ([Mni 2 Oi2(Ac)i 6 (H20) 4 ]-2HAc-4H20, where 
HAc=accetic acid), or Mni2 for short^. It usually has 
spin S = 10 and large anisotropy energy, producing 
a high spin reversal barrier—. Many interesting phe- 
nomena have been observed, such as various dynami- 
cal magnetism. One of the most intriguing phenom- 
ena observed in SMM systems is a step-wise structure 
in low-temperature magnetization curves^—. Great ef- 
forts have been made to investigate this phenomenon and 
related effects! 3 -^. The step -wise structure is attributed 
to Landau-Zener (LZ) quantum tunneling effec t 21 ' 22 . 
This stimulates intensive study on LZ model and its 
variants^ 3 - - — . Some authors use numeric diagonalization 
method a 32 : 33 to study many-level LZ models to under- 
stand the step structure in experimental magnetization 
curves. However, it is difficult to consider thermal effects 
in these approaches to obtain satisfactory magnetization 
curves comparable to experimental results. 

In this paper, we shall combine the classical thermal 
effects with the quantum LZ tunneling effects in a dy- 
namical Monte Carlo (DMC) framework 3 -^ - — in order 
to produce satisfactory magnetization curves compara- 



ble to experimental results. We consider ideal tetragonal 
body-centered lattices and use the giant spin approxima- 
tion for spins of SMMs. We consider magnetic dipolar 
interactions, but neglect other factors such as defects, 
disorders, and misalignments between the easy axis and 
applied magnetic field. We calculate spin reversal prob- 
abilities from thermal-activated barrier hurdling, direct 
LZ tunneling effect, and thermal-assisted LZ tunneling 
effects in the presence of sweeping magnetic fields, and 
thereby derive a unified probability expression for any 
temperature and any sweeping field. Taking the Mni2 as 
example, we do systematical DMC simulations with vari- 
ous temperatures and sweeping rates. The step structure 
appears in our simulated low-temperature magnetization 
curves, and our simulated magnetization curves are semi- 
quantitatively consistent with corresponding experimen- 
tal results on those good Mni2 systems (with less dis- 
orders) in the presence of little misalignments between 
the easy axis and applied field o 15 ' 16 . Interplays of the LZ 
tunneling effect, the thermal effects, and the magnetic 
dipolar interactions are elucidated. These imply that 
our simple model and DMC method capture the main 
features of experimental magnetization curves for little 
misalignments. More detailed results will be presented 
in the following. 

The rest of this paper is organized as follows. In next 
section we shall define our spin model and describe ap- 
proximation strategy. In the third section we shall de- 
scribe our simulation method, present our unified proba- 
bility formula for the spin reversal from the three spin re- 
versal mechanisms, and give our simulation parameters. 
In the fourth section we shall present our simulated mag- 
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netization curves and some analysis. In the fifth section 
we shall show the key roles of the dipolar interactions in 
determining LZ tunneling probabilities. Finally, we shall 
give our conclusion in the sixth section. 



II. SPIN MODEL AND APPROXIMATION 

Without losing generality, we take typical Mn 12 system 
as our sample in the following. Under giant spin approxi- 
mation, every Mn 12 SMM is represented by a spin 5=10. 
Magnetic dipolar interactions are the only inter-SMM in- 
teractions, with hyperfme interactions neglected. Mni2 
SMMs are arranged to form a body-centered tetragonal 
lattice with experimental lattice parameters 37 . Using a 
body-centered tetragonal unit cell that consists of two 
SMMs, we define our lattice as L\ x L 2 x L3, where L\, 
L2, and L3 are three positive integers. A longitudinal 
magnetic field B z {t) = Bq + vt is applied along the o 
easy axis of magnetization, where v is the field-sweeping 
rate and Bq is the starting magnetic field. The total 
Hamiltonian of this system can be expressed as 

where Hf is the single-body part for the i-th single SMM, 
and Hfj describes the magnetic dipolar interaction be- 
tween the i-th and j-th SMM. The factor 1/2 before the 
sum sign is due to the double counting in the summation. 
Hf is given by 

H?= -D{Stf+E[{Sff-{S y t f] 

+B° 4 6 4 + B\6\ + gf x B B z Sf , (2) 

where Si = (Sf , Sf , Sf) is the spin vector operator for the 
i-th SMM, g the Lande g- factor (here g — 2 is used), [i B 
the Bohr magneton, D, E, B\ and B\ are all anisotropic 
parameters, and 0\ and 0\ are both Steven operators^ 
defined by 0° = 35(5? ) 4 - [305(5 + 1) - 25](5?) 2 + 
35 2 (5 + l) 2 - 65(5 + 1) and 6% = [(5+) 4 + (5r) 4 ]/2. 
Hf- is defined by 

Hfl = ^#[^ • S, - -^(S, • • m)] , (3) 

ij ij 

where /zo is the magnetic permeability of vacuum, and r,j 
the vector from i to j, with rjj=|ry| being the distance 
between i and j. 

For the i-th SMM, we treat all the effects from the 
other SMMs by classical-spin approximation. As a result, 
we derive the partial Hamiltonian Hi that acts on the i-th 
SMM: 

Hi = H? + g^ B Bf ■ S, 

= -D{Stf + BlOl + Hf + g^ B {B z + B%)Sf ,(4) 



where the transverse part Hf is defined as 

Hf = E[{Sff-(Sm + Bt6j+gMB%S* + B%S!).(5) 

For the i-th SMM, the dipolar interaction of the 
other SMMs is equivalent to Bf = {Bg,B^,Bg) = 
J2j(=ii) Bjjj where is the magnetic dipolar field ap- 
plied by the j-th SMM on the i-th SMM. It contributes 
a magnetic field consisting of longitudinal and transverse 
parts. 

III. SIMULATION METHOD AND 
PARAMETERS 

As we show in Fig. 1, there are three main mechanisms 
related to the reversal of a SMM spinST. 1 ^^^: ( a ) 
thermal-activated barrier-hurdling, (b) direct LZ tunnel- 
ing, and (c) thermal-assisted LZ tunneling. The thermal- 
activated barrier hurdling dominates at high temperature 
(if the blocking temperature Tb ^3.3K for Mni^ is 
treated as high temperature), the direct LZ tunneling at 
low temperature, and the thermal- assisted LZ tunneling 
at intermediate temperature. For any temperature, we 
consider all the three spin reversal mechanisms simulta- 
neously. For the time scale we are interested, we do not 
need to treat phonon-related interactions directly, but 
shall use an effective transition-state theory to calculate 
the thermal-activated spin-reversal rates. We shall use 
a DMC method to combine the quantum LZ tunneling 
effects with the classical thermal effects. Various kinetic 
Monte Carlo (KMC) methods^—, essentially similar to 
this DMC method, have been used to simulate atomic 
kinetics during epitaxial growth for many years. On 
the other hand, MC simulation has been used to study 
Glauber dynamics of kinetic Ising models^—. We shall 
present a detailed description of this DMC simulation 
method in the following. 

A. Thermal-activated spin reversal probability 

We need the thermal-activated energy barrier in or- 
der to calculate the thermal-activated spin-reversal rate. 
When calculating the thermal-activated energy barrier 
we ignore the small transverse part Hf and use classi- 
cal approximation for the spin operators. The large spin 
5 = 10 of Mni2 further supports the approximations. As 
a result, the energy of the i-th SMM can be expressed as 

Ei = -D 2 {Stf - D 4 {S*f + hiS? , (6) 

where Sf is the classical variable for the spin operator 
Sf, hi = gf i B {B z + Sg), D 2 = D + [305(5+1) - 25]S|, 
and D4 — — 35-B2. Because hi is dependent on time t, Ei 
changes with t. 

We define our MC steps by the time points, t n = At ■ n, 
where n takes nonnegative integers in sequence. For the 
n-th MC step, we use Ei >n , hi^ n , and Sf n to replace Ei, 
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FIG. 1: (Color online.) A schematic demonstration of the three spin reversal mechanisms: (a) thermal-activated barrier 
hurdling, (b) direct LZ tunneling, and (c) thermal-assisted LZ tunneling. The probabilities, energy levels, barrier, and other 
symbols are defined in the text. The horizontal solid line with arrow in (b) and (c) shows that a pair of energy levels satisfy 
the resonance tunneling conditions. The horizontal dotted lines in (b) and (c), as guide for eyes, imply that these energy levels 
do not match. 



hi, and Because each of the spins has two equi- 

librium orientations along the easy axis, we assume ev- 
ery spin takes either S or — S at each of the times t n . 
Within the n-th MC step (t: from t n to t n+ i), we use 
an angle variable 0^ n to describe the i-th spin's devi- 
ation from its original (t n ) orientation S*^. Naturally, 
Qi,n — corresponds to the original state and Qi iU — tt 
is the reversed state, and then all the other angle values 
(0 < 0i <n < 7t) are treated as transition states. Express- 
ing Sf n as S*^ cos 9i tn , we usually have a maximum in 
the curve of E^ n (cos 9i jn ) as a function of cos 6^ n , and the 
maximum determines the energy barrier for the spin re- 
versal mechanism^— , as shown in Fig. 1(a). We define 
^j,n = cos0i in for convenience. We have — 1 < xi iTl < 1 
for actual 9^ n , but be extended beyond this re- 

gion in order to always obtain a formal solution x™^ x for 
the maximum. |xf^ x | < 1 implies that there actually ex- 
ists an energy barrier, and |a;f^ x | > 1 means that there 
is no barrier for the corresponding process. Under con- 
ditions D-2 > and D4 > 0, the barrier can be expressed 
as: 

r ^,n«r)> Kn x \ < 1 

AE i>n = lE i>n (-l) = \2hi, n S^ n \, <T<-1 (7) 
{ E i<n (l) = 0, > 1 

where xf 1 ^ is defined by 

<n x = y/~qi, n /2+y/d~+ y/-qi, n /2 - yfd~ (8) 

and the three parameters are defined by dj „ = „/2) 2 + 
(p/3) 3 , p = D 2 /(2D 4 S 2 ), and q i<n = -h^S^J (4D 4 S 4 ). 
These parameters are dependent on the spin configura- 
tion and the magnetic field, and then on the time t n (or 
n). 

The spin reversal rate within the n-th MC step 
(between t n and t n +i) can be expressed as Ri^ n = 
Ro exp(— AEi^n/ksT) in terms of Arrhenius law^, where 



kg is the Boltzmann constant and R the characteristic 
attempt frequency. We use P n (t') to describe the prob- 
ability that the i-th spin is reversed between and t' , 
where t' satisfies the condition t' < At. It has the ini- 
tial condition P n {t' = 0) = and satisfies the equation 
[1 - P n {t')] ■ R n (t')dt' = P n (t' + dt') - P n (t'), or 

[1 - P n (t')]R n (t') = ^Pn{t>) (9) 

where R n (t'), the reversal rate at t' , is taken as the rate 
Ri,m independent of t' within the region [0,Ai]. Solving 
the equation, we obtain the probability Pf^ defined as 
P n (t' = At) for a classical thermal-activated reversal of 
the i-th spin within the n-th MC step: 

P^ s = l-exp(-At-R^ n ). (10) 

For At < 1/Ri,n, Eq. (JDJl reduces to = At ■ FU >n . 
The probability expression defined in Eq. (fTUj) is reason- 
able because Pfjf 55 will not exceed unity even when At is 
very large with respect to l/i?i iK . 



B. LZ-tunneling related spin reversal probabilities 

When temperature is lower than Tb , LZ tunneling be- 
gins to contribute to spin reversal. We begin with the 
effective quantum single-spin Hamiltonian (jj) with (JSJ. 
All the effects of other spins are included in the mag- 
netic dipolar field Bf 1 (depending on the time t) and 
are depending on the magnetic field and the current spin 
configuration. For the n-th MC step, if the transverse 
term H\ T is removed, Hamiltonian Eq. (Q} is diagonal 
and has 25+1 energy levels, E^, where m can take 
any of S, S — 1, • ■ • , —(S — 1), — S. If using the continuous 
time variable t, we can express the energy levels as E l m (t) 
(with m from S to -S) and derive their crossing fields [at 
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which El{t) =E\ n , (t)] : 

(m + m')[L» 2 + £>4(to 2 + m' 2 )] 



P> 



5MB 



(11) 



The transverse term Hf T will modify the energy levels 
fijf , but the 25 + 1 energy levels of Hamiltonian Q with 
©, P4 n , can be still labelled by m = S, 5-1, • • • , -(5- 
1), — 5. Actually, the difference between E]£ and P^ ra 
is small. Due to the existence of the transverse part 
P* r , there will be an avoided level crossing between P^™ 
and E 1 ^ for the n-th MC step when the effective field 
B z + Bfl equals B % ^ m ,, with m and m' taking values 

among 5, 5-1, • • • , -(5-1), -5. The set of all the B^ m , 
values are the effective field conditions for the avoided- 
level-crossings. If E 1 ^ 1 equals P^?, B % ^ m , is approxi- 
mately equivalent to the crossing field (equaling P m , m '). 
The allowed (m, m') pairs are shown in Fig. 2. This 
means that when B z is swept to a right P^™ m / — Bf* 
value, a quantum tunneling occurs between the m and 
to' states. The tunneling can be well described using 
LZ tuimelingi£i22i22i2£. The nonadiabatic LZ tunneling 



probability P 



LZ, 



is given bya»: 



21,22 



P 



LZ,j,i 



1 — exp 



tt(a: 



».n \2 
m,m' / 



2hg/iB\m — m!\v 



(12) 



where the tunnel splitting m , is the energy gap at the 
avoided crossing of states m and to'. B % 2^ , and A!^ m / 
can be calculated by diagonalizing Eq. (J3J ■ If the dipo- 
lar field is neglected, P^" m „ A^ m „ and P^;" reduce 
to P^. m /, A^ m ,, and P^ tm ,, those of corresponding iso- 
lated SMMs, respectively. 

At the beginning of field sweeping, we let all the spins 
have to = 5. If T <^Tb, thermal activations are frozen, 
and LZ tunnelings only occur at the avoided crossings 
(5, to'), where to' takes one of —5, —5 + 1, • • • , 5 — 1. 
This is the direct tunneling shown in Fig. 1(b), and the 
LZ tunneling probability is given by 



pdLZ pdir,i. 

i,n S,m' 



(13) 



It is nonzero only when the condition Ef 1 — E^J is 
satisfied. When the temperature is in the intermediate 
region <C T < Tb , the thermal-assisted tunneling plays 
an important role. This process can be represented by 
5 to — > m' as shown in Fig. 1(c), in which 5 and 
to states lie on one side of the thermal barrier and to' 
and —5 states on the other side. The first process 5 
to means that a spin is thermally activated from 5 to 
to state with the probability Pg^'^™, which is given by 
p S C ^m n = 1 - exp(-At • Rf%), where Pf^ is given by 
Poexp[— (P^™ — Eg n )/kBT]. The second process m — > 
in' is the LZ tunneling from m to m! ', with the probability 
defined in Eq. (TT2"|) . Therefore, the reversal probability 



of thermal-assisted LZ tunneling through m is given by 

(14) 



ptaLZ pass,i,n pact,i,n nLZ,i,n 

*i,n,m S-~>m—¥m' S-~>m m,m' ' 



It is nonzero only when the condition P^™ 
satisfied. 

It must be pointed out that to' in P s 



and 



, is determined by P^™ 



P^andP^P^, 



respectively, as is shown in Figs. 1(b) and 1(c). If the 
energy-level condition is satisfied, the probability is larger 
than zero; or else the probability is equivalent to zero. 
Therefore, the subscript to' in Pg 1 ^!'™ and P^l^ m , can 
be removed, as we have done in P^ z and Pt^Jk- 



Unified spin reversal probability for MC 
simulation 



Generally speaking, every one of the three spin rever- 
sal mechanisms takes action at any given temperature. 
Actually the LZ tunneling effect dominates at low tem- 
peratures and the thermal effects become more important 
at higher temperatures. For the n-th MC step, the prob- 
ability for the thermal-activated barrier-hurdling reversal 
of the i-th spin is given by P^ s defined in Eq. ijTUjl [sec 
Fig. 1(a)], that for the direct LZ tunneling effect equals 
P^ z defined in Eq. (USD [see Fig. 1(b)], and that for 
the thermal-assisted LZ tunneling effects through the m 
state is given by P*^ defined in Eq. (0 [see Fig. 1(c)]. 
Here the partial probabilities from the three mechanisms 
are considered independent of each other. Therefore, we 
can derive the total probability P*°' for the reversal of 
the i-th spin within the n-th MC step: 

p^=i-(i-p^)(i-p^ z ) n 

mtop <m<S 

(15) 

where m top , depending on the effective field, is deter- 
mined by the highest level among the 25 energy 

levels, P^™ (-5 < to < 5), as we show in Fig. 1(c). 

It must be pointed out that Pf„ is always larger than 
zero, but the LZ-tunneling related probabilities, P^ z 
and P'^ L m, are nonzero only at some special values of 
the effective field. As is shown in Fig. 2, there is at most 
one LZ-tunneling channel, from either direct or thermal- 
assisted LZ effect, for a given nonzero value of the effec- 
tive field. As a result, when the effective field is nonzero, 
we have at most one nonzero value from either Pf^ or 
one of P^nm ( m to P < m < 5). It is only at the zero 
value of the effective field that both P^ z and P*^ LZ 
(0 < m < 5) (m top = 0) can be larger than zero so that 
we can have the direct LZ tunneling and all the thermal- 
assisted LZ-tunneling channels simultaneously. In our 
simulations, the processes that a reversed spin is reversed 
again are also considered, but the probabilities are tiny. 
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FIG. 2: (Color online.) A schematic demonstration for the conditions of the i-th spin that Landau-Zener tunnelings can happen. 
m labels the spin z component, from -10 to 10, and B z + Bfl is the effective magnetic field. A hollow circle indicates one of the 
allowed m values. The two circles (m and m') connected with one vertical line means that a LZ tunneling condition is satisfied 
between the two states at the corresponding effective field S*!" _, within the n-th MC step. There is at most one LZ tunneling 
at any nonzero value of the effective field, but every m state can tunnel to the -m state when the effective field is zero. The 
inset amplifies the part between 0.3 and 1.2 T. 



D. Simulation parameters 

We use experimental lattice constants, a = b = 
17.1668 A and c = 12.2545 A, and experimental 
anisotropy parameters, D/ks = 0.66K, B^/ks = 
-3.2 x 10~ 5 K, and Bf/k B = 6 x 10~ 5 KiSiiML. As 
for the second-order transverse parameter E, Ejks = 
1.8 x 10 _3 K is taken from the average of experimental 
values^. We describe the time by using both continuous 
variable t and discrete superscript /subscript n. In some 
cases, the sweeping field can be used to describe the time 
because it is defined by B z (t) = Bo + ft. There is al- 
ways a nonnegative integer n for any given t value, and 
there is a t region, [t n ,t n+ i], for any given nonnegative 
n. We take At — 0.1ms and i?o = 10 9 /s, which guar- 
antee the good balance between computational demand 
and precision. 

The dipolar fields (B%, B$, Bg) at each SMM are up- 
dated whenever any of the SMM spins is reversed. The 
^Knm' va l ues are recalculated whenever any LZ tunnel- 
ing happens. In the simulations, the field B z is swept 
from -7 to 7 T in the forward process, and the full mag- 
netization hysteresis loop is obtained simply by using 



the loop symmetry. Every magnetization curve is cal- 
culated by averaging over 100 runs to make statistical 
errors small enough. The main results presented in the 
following are simulated and calculated with lattices con- 
sisting approximately of 900 ~ 1200 body-centered unit 
cells or 1800 ~ 2400 spins. We have test our results with 
lattices consisting approximately of 100 ~ 6000 body- 
centered unit cells, or 200 ~ 12000 spins. 



IV. SIMULATED MAGNETIZATION CURVES 

Presented in Fig. 3 are simulated magnetization curves 
(with M normalized to the saturated value Ms) against 
the applied sweeping field B z for ten different tempera- 
tures: 0.1, 0.5, 0.6, 0.8, 1.0, 1.5, 2.0, 2.5, 2.8, and 3.2 K. 
Here, the lattice dimension is 10 x 10 x 10 and the field 
sweeping rate is 0.02 T/s. Each of the curves is calculated 
by averaging over 100 runs. The curves of 0.1 K and 0.5 K 
fall in the same curve, which implies that thermal acti- 
vation is totally frozen when the temperature is below 
0.5 K. It can be seen in Fig. 3 that the area enclosed 
by a magnetization loop decreases with the temperature 
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increasing, becoming nearly zero at 3.2 K (near the block- 
ing temperature 3.3 K of Mn 12 ). There are clear magne- 
tization steps when the temperature is below 2.0 K. They 
are caused by the LZ quantum tunneling effects. For con- 
venience, we describe a step by using a H-part, a vertex, 
and a V-part. For an ideal step, the H-part is horizon- 
tal and the V-part vertical, but for any actual step in 
a magnetization curve, the H-part is not horizontal and 
the V-part not vertical because of the dipolar interaction 
and thermal effects, and the two parts still meet at the 
vertex. The vertex is convex toward the up-left direction 
in the right part of a magnetization loop and toward the 
down- right direction in the left part. At higher temper- 
atures (> 2. OK), there is no complete step and there are 
only some kinks that remind us of some LZ tunnelings. 
This should be caused mainly by thermal effects. 



LZ tunneling effects. Different sweeping rates lead to 
substantial changes in the magnetization curves, and the 
larger the sweeping rate becomes, the larger the hystere- 
sis loops are. 




FIG. 3: Simulated magnetic hysteresis loops (M/Ms vs B z ) 
with sweeping rate v — 0.02 T/s for ten temperatures: 0.1, 
0.5, 0.6, 0.8, 1.0, 1.5, 2.0, 2.5, 2.8, and 3.2 K (from outside to 
inside). The lattice dimension is 10 x 10 x 10. Note that the 
two curves of 0.1 K and 0.5 K fall in the same curve. 

Presented in Fig. 4 are the right parts of the mag- 
netization curves against the applied sweeping field for 
three temperatures, 0.1, 1.5, and 2.5 K, and with three 
sweeping rates, 0.002, 0.02, and 0.2 T/s. Here, the lat- 
tice dimension is 10 x 10 x 10. We label a magnetiza- 
tion step by the magnetic field defined by its V-part near 
its vertex. For T=0.1K, only the direct LZ tunnelings 
change the magnetization, and the magnetization steps 
from B z =2 to 6T in Fig. 4 correspond to B° s , with m' 
being from -6 to 2 in Table I. For T=1.5K, there are clear 
steps in the lower parts of the three magnetization curves, 
but their V-parts deviate substantially from the corre- 
sponding Bg , values and the steps are substantially 
deformed, which show that thermal-assisted LZ tunnel- 
ings play an important role. When temperature rises to 
2.5K, there is no step structure and only one kink can 
be seen in the lower part of the magnetization curve in 
the cases of 0.2T/s and 0.02T/s. This is because the ef- 
fects of thermal activation become dominating over the 




FIG. 4: (Color online.) The right parts of simulated mag- 
netization curves of different sweeping rates 0.002 (dot), 0.02 
(dash), and 0.2 (solid) T/s for three temperatures 0.1 K, 1.5 K, 
and 2.5 K, as labelled. The lattice dimension is 10 x 10 x 10. 
Each of the visible steps and kinks along a magnetization 
curve corresponds to one of the magnetic fields at which the 
direct and thermal assisted LZ tunnelings take place. The 
thin vertical dotted lines show the positions of B° s m , for ra —- 
10,-9,- -.,2. 




B z (T) 

FIG. 5: (Color online.) The right parts of simulated mag- 
netization curves for three temperatures with five different 
lattice dimensions: 20 x 20 x 3 (dash-dot), 12 x 12 x 8 (dash), 
10x10x10 (solid) ,9x9x14 (dot) , and 3 x 3 x 100 (short-dash) . 
The temperatures are 0.1 K, 1.5 K, and 2.5 K, as labelled. The 
sweeping rate is 0.02 T/s. For comparison, we also present 
the results without considering dipolar interaction (thin solid 
line). 

Presented in Fig. 5 are the right parts of simulated 
hysteresis loops with i/=0.02 T/s at three temperatures 
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for five different lattice dimensions: 20x20x3, 12x12x8, 
10 x 10 x 10, 9 x 9 x 14, and 3 x 3 x 100. The temperatures 
are 0.1, 1.5, and 2.5 K. For comparison, the simulated re- 
sults without the dipolar interaction are presented too. 
For T=0.1K, there are clear step structures for all the 
five lattice shapes. The step height varies with the lattice 
shape, which can be attributed to the dipolar interaction. 
If the dipolar interaction is switched off, there are only 
two steps: one tall step at B z — 4.00 T and the other very 
short step at 3.06 T. They correspond to the two transi- 
tions from 10 to -2 and -4, respectively. Other transitions 
from 10 to -4, -6, -8, and -10 have too small probabilities 
to be seen. When the dipolar interaction is switched on, 
the transition from 10 to -3 is allowed and the tall step be- 
comes much shorter, resulting in the rich step structures 
between 3 and 6 T. The steps are caused by the direct 
LZ tunnelings. When the temperature changes to 1.5K, 
the hysteresis loops become substantially smaller because 
of the enhanced thermal effects. In this case, there are 
deformed step structures in the lower parts of the mag- 
netization curves and there does not exists any clear step 
structure in the upper parts. The deformed step struc- 
tures between 1 and 3 T result from the thermally as- 
sisted LZ tunnelings. For T=2.5K, there does not exist 
any step structure at all for all the six cases. The effect of 
the lattice shape is attributed to the long-range property 
of the dipolar interaction, and can be clearly seen in the 
magnetization curves only at the low temperatures in the 
extreme cases of 20 x 20 x 3 and 3 x 3 x 100. Actually, 
there is little visible difference between the magnetiza- 
tion curves of the three lattices: 12 x 12 x 8, 10 x 10 x 10, 
and 9 x 9 x 14. Visible difference can be found at 0.1K 
and 1.5K only for the two extreme cases: 20 x 20 x 3 and 
3 x 3 x 100. If we define a ratio r = Li/L t of longitudinal 
size to transverse size for L t x L t x L;, we have r=l for 
10 x 10 x 10, r=0.67 for 12 x 12 x 8, r=1.56 for 9 x 9 x 14, 
r=0.15 for 20 x 20 x 3, and r=33 for 3 x 3 x 100. There- 
fore, there is little clear effect of lattice shape as long 
as the shape parameter r is neither extremely large nor 
extremely small. 

Now we address the statistical errors. We have cal- 
culated standard errors <tm of the reduced magnetiza- 
tion M/M s as functions of the sweeping field for various 
temperatures and sweeping rates. Our results show that 
for a given magnetization curve, the statistical errors are 
very small (<7m < 0.005) in the region of B z defined by 
\M/M S \ > 0.9, and reach a maximal value cr^ ax near the 
point of B z defined by M/M s — 0. The maximal sta- 
tistical error is dependent on the temperature and 
sweeping rate, varying from 0.015 to 0.025 for our simu- 
lation parameters. Such statistical errors appear only in 
a very small region of B z . For any magnetization curve 
as a whole, the statistical errors are small enough to be 
acceptable. 

Here we discuss effects of lattice sizes on simulated 
results. The above simulated results are based on the 
lattice dimensions: 20 x 20 x 3, 12 x 12 x 8, 10 x 10 x 10, 
9 x 9 x 14, and 3 x 3 x 100. They have 900 - 1200 body- 



centered unit cells, or 1800 ~ 2400 spins. To test our re- 
sults, we have done a series of simulations for different pa- 
rameters using lattice dimension defined by L t x L t x Li. 
In the cases of T = 0.1K, v = 0.2T/s, and L t = L t = L 
with L = 5 ~ 20, the largest size effects appear between 4 
and 5.5 T for the right parts of the magnetization curves. 
For the steps at 4T, the L-caused change in the magne- 
tization decreases quickly with increasing L, becoming 
very small when L is larger than 9. Therefore, our lat- 
tice sizes of the results presented above are large enough 
to be reliable. 

The above simulated results show that the area en- 
closed by a magnetization hysteresis loop decreases with 
the temperature increasing and increases with the sweep- 
ing rate increasing. This is completely consistent with 
the temperature and sweeping-rate dependence of the 
thermal reversal probability and LZ tunneling probabil- 
ities. Thermal activation effects dominate at high tem- 
perature. The LZ tunneling effects manifest themselves 
through the steps and kinks along the magnetization 
curves. However, there is a limit for the hysteresis loops 
at the low temperature end for a given sweeping rate. 
These limiting magnetization curves are caused by the 
minimal reversal probability set by the direct LZ quan- 
tum tunneling effect because the thermal activation prob- 
ability becomes tiny at such low temperatures. With 
usual shape parameter r, these results are consistent with 
experimental magnetization curves of good Mni2 crystal 
samples in the presence of little misalignments between 
the easy axis and applied field o 15 ' 16 . In principle, a trans- 
verse magnetic field (due to the misalignment of the ap- 
plied field and the easy axis) can enhance the energy 
splitting, and as a result will reduce the magnetization 
loop and smooth some step a 15 ' 16 ' 50 . These usual (not 
extreme) shape parameters should reflect real shape fac- 
tors in experimental samples. The consistence should be 
satisfactory, especially considering that our theoretical 
probabilities are calculated under leading order approxi- 
mation and our model does not include possible defects 
and disorders in actual materials. 



V. KEY ROLES OF DIPOLAR FIELDS 

To investigate the effects of dipolar interactions, we 
divide the dipolar fields within the n-th MC step, 
(Bf^ n , Bfy n , Bf^ n ), into two parts: transverse dipolar 
field Bf* n and Bf* n , and longitudinal dipolar field Bfl n . 
Transverse dipolar field not only modifies B^f m ,, but 
also affects AY 1 , and -P^ 2 ^'". In contrast, longitudinal 
dipolar field affects neither A^" m , nor P^ 2 ^'", but shifts 
-B™"™/ by —Bfl „ . This means that LZ tunnelings actu- 

ally occur at the field B^ m , - B?> nl not B l ^ m ,. This 
shift has two effects. First, it broadens the LZ transi- 
tion and deforms the steps in magnetization curves. Sec- 
ond, the quick changing of Bfl results in that the value 



TABLE I: Calculated results of B° Sm ,, 



p o 



{P Sr ff}, and <Tg m , for the direct LZ tunneling (S,m') when 
the field B z is swept to 3.75 T, where n is determined by 
the field 3.75 T. T = 0.1 K, v = 0.02 T/s, and the lattice 
dimension is 10 x 10 x 10. 



A£g, m ,(T) 



p o 



-10 
-9 
-8 
-7 
-6 
-5 
-4 
-3 
-2 
-1 

1 
2 
3 



0.000000 
0.564160 
1.099966 
1.612415 
2.106511 
2.587260 
3.059671 
3.528757 
3.999529 
4.476997 
4.966165 
5.472035 
5.999604 
6.553867 



6.4xl0" 15 

1.6xl0~ 6 

3.5xl0" 6 

5.1xl0~ 6 

6.7xl0~ 6 

7.9xl0 -6 

8.6xl0~ 6 

8.6xl0 -6 

7.8xl0 -6 

6.3xl0~ 6 

3.9xl0~ 6 

7.4xl0~ 7 

3.6xl0~ 6 

8.6xl0~ 6 



/ pLZ,n\ 



' S,m' 



0.00000 


0.00000 


0.00138 


0.01815 


1.00000 


1.00000 


1.00000 




0.00000 
0.00000 
0.00000 
0.00000 
0.00138 
0.00002 
0.01838 
0.22194 
1.00000 
0.53746 
1.00000 
0.99975 
1.00000 
0.99988 



0.00000 
0.00000 
0.00000 
0.00000 
0.00001 
0.00002 
0.00320 
0.21086 
0.00000 
0.33455 
0.00089 
0.01091 
0.00000 
0.00749 



Blnm' can b e m i sse d by the effective field B z + Bf* and 
therefore the actual percentage of the reversed spins due 
to the LZ tunneling effect with respect to the total spins 
is smaller than the LZ probability -P™ given in Eq. 
(fl"2|) . This means that the dipolar interaction hinders 
both the direct LZ tunneling process and the thermal 
assisted LZ tunneling processes. 



Without transverse dipolar field, B % ^ m , becomes 

Bm m< j an d Ps'm'' n ecma ls for odd m! values because 
transverse dipolar field is the only transverse term of 
odd order in Hamiltonian Eq. (4). Without longitu- 
dinal dipolar field, the V-parts of steps remain vertical 
and the percentage of the reversed spins due to LZ tun- 
neling is strictly equivalent to the LZ probability Pg J^i' n 
at low temperatures. These arc shown by the thin solid 
line for 0.1K in Fig. 5. In Table I we also present the av- 
erage value (AB% tm ,=(\B§ tm , ~ Bg iTB ,|» of dipolar-field 
fluctuations with respect to B Sm ,, the dipolar- field- free 
LZ probability Pg m >, and the average value (Pg Z '' 



S,m> I aild 
LZ,i,n 
S,m' Ui r S,m> 

avoided crossing positions of S and to', where to' varies 
from -10 to 3 and the averaging (X n ) of X l,n is calcu- 
lated over all the spins and all the runs within the n-th 
MC step. It should be pointed out that the ABg , val- 
ues, although very important to LZ tunnelings, are very 
small, as shown in Table I. It is transverse dipolar field 
that make ABg m , nonzero and make Pg^ 1 /" (m'=-5, -3, 
-1, 1, 3) change from to nonzero, even nearly reach 1 
in the cases to' = 1 and 3. 

In order to elucidate the magnitude and distribution of 
the dipolar fields, we address the time-dependent distri- 



the corresponding standard error <Jg m , of P s ^!,' n for the 



75000 
60000 
45000 
30000 
15000 


30000 
20000 
10000 


30000 
20000 
10000 


30000 
20000 
10000 


100000 
75000 
50000 
25000 




1 (a) 20x20x3 / 

/ 
/ 

- j$ 


P ' di - 

\ B ix ■ 

k ^ag — : 


- (b) 12x12x8 $ 


a i i i i - 

^ : 
»V - 

I 1 ■ • . i gma.fflna.aA 


1 1 1 1 
- (c) 10x10x10 $ 


1 1 1 1 

Jg - 


1 1 1 1 5 
- (d) 9x9x14 

j y g. Q □ ^ q □ ■ 1 


s 1 1 1 1 - 
h - 


- 1 1 1 1 fl 

■ (e) 3x3x100 

j 

A / 


1 




\ 



-0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 0.1 
Dipolar field (T) 

FIG. 6: (Color online.) Distributions of dipolar fields Bf^ 
(dashed line + circle), Bfy (dotted line + cross), and Bfl 
(solid line + square) for the five lattice dimensions 20 x 20 x 3 
(a), 12 x 12 x 8 (b), 10 x 10 x 10 (c), 9 x 9 x 14 (e), and 
3 x 3 x 100 (e) when the field B z is swept to 3.75 T. The 
temperature T is 0.1 K and the sweeping rate v equals 0.02 
T/s. 



butions of SMMs that have dipolar fields {Bf>, Bf>, Bf>) 
(here the continuous time variable is implied), or in short 



the distributions of Bf*, Bf^, 



and Bf*, in the following. 



In Fig. 6 we compare the results from five different lat- 
tice dimensions: 20 x 20 x 3, 12 x 12 x 8, 10 x 10 x 10, 
9 x 9 x 14, and 3 x 3 x 100. Here the time is when the 
field B z is swept to 3.75 T, the temperature T is 0.1 K, 
and the sweeping rate v equals 0.02 T/s. For all the five 
lattices, our results show that the distribution of Bf* al- 
ways is approximately equivalent to that of Bfy and they 
are both symmetrical and peaked at zero. The peak is 
sharper for the extremely slab-like 20 x 20 x 3 lattice 
and extremely rod-like 3 x 3 x 100 lattice. The peak 
of the Bfl distribution is wider than that of both Bf^ 
and Bfy. It shifts substantially away from zero when the 
lattice shape is either extremely slab-like or extremely 
rod-like. The leftward shift of the Bfl peak can be at- 
tributed to dipolar-interaction-induced ferromagnetic or- 
ders in rod-like system a 51 ' 52 , and the similar rightward 
shift to dipolar-interaction-induced antiferromagnetic or- 
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ders in slab-like systems. Because dipolar interactions 
are the only inter-SMM interactions in our model, the 
differences of distributions between the the five lattices 
are caused by the dipolar fields, or dipolar interactions 
in essence. 

VI. CONCLUSION 

In summary, we have combined the thermal effects 
with the LZ quantum tunneling effects in a DMC frame- 
work by using the giant spin approximation for spins of 
SMMs and considering magnetic dipolar interactions for 
comparison with experimental results. We consider ideal 
lattices of SMMs consistent with experimental ones and 
assume that there are no defects and axis-misalignments 
therein. We calculate spin reversal probabilities from 
thermal-activated barrier hurdling, direct LZ tunneling 
effect, and thermal- assisted LZ tunneling effects in the 
presence of sweeping magnetic fields. Taking the param- 
eters of experimental Mni2 crystals, we do systematical 
DMC simulations with various temperatures and sweep- 
ing rates. Our results show that the step structures 
can be clearly seen in the low-temperature magnetiza- 
tion curves, the thermally activated barrier hurdling be- 
comes dominating at high temperature near 3K, and the 
thermal- assisted tunneling effects play important roles at 



the intermediate temperature. These are consistent with 
corresponding experimental results on good Mni2 sam- 
ples (with less disorders) in the presence of little mis- 
alignments between the easy axis and applied field o 15 ' 16 , 
and therefore our magnetization curves are satisfactory. 

Furthermore, our DMC results show that the magnetic 
dipolar interactions, with the thermal effects, have im- 
portant effects on the LZ magnetization tunneling effects. 
Their longitudinal parts can partially break the reso- 
nance conditions of the LZ tunnelings and their trans- 
verse parts can modify the tunneling probabilities. They 
can clearly manifest themselves when the SMM crystal is 
extremely rod-like or slab-like. However, both the mag- 
netic dipolar interactions and the LZ tunneling effects 
have little effects on the magnetization curves when the 
temperature is near 3K. This DMC approach can be ap- 
plicable to other SMM systems, and could be used to 
study other properties of SMM systems. 
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